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摘要 : 物种 的 遗传 多 样 性 是 物种 适应 性 和 生存 能 力 的 关键 因素 。 和 生境 片段 化 是 造成 生物 多 样 
性 丧失 的 重要 因素 之 一 , 其 对 植物 种 群 的 遗传 多 样 性 有 着 重要 影响 ,水 角 (Hydrocera triflora) 
处 于 地 区 性 绝 灭 状态 ， 其 遗传 多 样 性 状况 和 濒危 机 制 尚未 有 报到 。 该 文 收 集 了 7 个 种 群 的 
34 个 水 和 角 样 本 ， 利 用 简化 基因 组 测序 技术 (RAD-seq) 获得 了 单 核 背 酸 变异 位 点 (SNP) 。 
通过 种 群 遗 传 多 样 性 和 遗传 结构 的 分 析 , 并 结合 种 群 历史 动态 分 析 和 不 同 气候 情景 下 物种 潜 
在 分 布 区 预测 ， 探 讨 了 水 角 的 濒危 机 制 。 结 果 表 明 : (1) 水 角 遗 传 多 样 性 较 低 (Ho=0.1569、 
He-0.1654, Pi=0.1865) ， 遗 传 分 化 系数 较 高 ;， AMOVA 分 析 表 明 ， 遗 传 变异 主要 发 生 种 群 
内 ; (2) Mantel 检测 表明 环境 距离 与 遗传 距离 、 地 理 距 离 均 具 有 显著 的 正 相 关 关 系 ， 分 别 
为 P-0.044 7 和 P=0.008 2; (3) 水 角 的 有 效 种 群 大 小 从 全 新 世 中 期 开始 持续 下 降 ， 琼 北 火 
山 群 的 爆发 和 人 来 活动 起 着 不 可 忽视 的 影响 ， 〈4) 生态 位 模拟 表明 与 当代 气候 相 比 ， 在 未 
来 气候 变化 下 水 角 的 潜在 分 布 区 总 面积 变动 不 大 ， 但 在 高 CO;s 排 放 的 情景 下 ， 大 量 的 高 适 
生 区 将 会 丧失 并 转化 为 低 适 生 区 , 特别 是 位 于 马 来 群 岛 的 适 生 区 几乎 完全 消失 。 结果 表明 生 
境 片 段 化 导致 了 水 角 较 低 的 遗传 多 样 性 和 有 效 种 群 大 小 的 持续 下 降 。 因此， 自身 更 新 能 力 低 
以 及 人 为 活动 干扰 、 城 市 化 等 不 利 的 环境 条 件 是 导致 其 濒危 的 主要 原因 。 建议 加 强 对 水 角 的 
就 地 保护 ， 进 行人 工 授 粉 等 方法 提高 其 基因 流 以 增加 其 种 群 的 遗传 多 样 性 。 另 一 方面 ， 要 着 
重 保 护 湿 地 免 遭 破坏 。 
关键 词 : 生境 片段 化 ， 简 化 基因 组 测序 ， 遗 传 多 样 性 ， 遗 传 结构 ， 物 种 分 布 模型 ， 种 群 历史 
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Abstract: The genetic diversity of species was a key factor in their adaptability and survival 
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loss, and it had significant impacts on the genetic diversity of plant populations. Hydrocera 
triflora was in a regional extinction state, and the status of its genetic diversity and endangered 
mechanisms had not been reported. In this study, 34 samples of H. triflora from 7 populations 
were used to obtain single nucleotide variation sites by restriction site-associated DNA sequencing 
(RAD-seq). We discussed the mechanisms of endangerment of H. triflora by analyzing population 
genetic diversity and structure. Additionally, we combined these findings with historical 
population dynamics analysis and predictions of potential species distribution in different climate 
scenarios. The results were as follows: (1) H. triflora had low genetic diversity (Ho=0.1569, 
He=0.1654, Pi=0.1865) and low genetic differentiation; AMOVA analysis showed that genetic 
variation mainly occurred within populations; (2) The Mantel test indicated that there was a 
significant positive correlation between environmental distance and genetic distance, as well as 
geographical distance, with respective p-values of 0.0447 and 0.0082; (3) The effective population 
size of H. triflora had been continuously declining since the mid-Holocene. The eruption of the 
Qiongbei volcanic group and human activities played a significant role; (4) Compared to the 
modern climate, the total potential distribution area of H. triflora changed little in the future. 
However, in scenarios of high CO» emissions, the highly suitable area decreased significantly, 
especially in the Malay Islands where the suitable habitat virtually disappeared. The results 
indicated that habitat fragmentation led to a continuous decline in the genetic diversity and 
effective population size of the H. triflora. Therefore, the low self-renewal capacity, as well as 
detrimental environmental conditions such as human disturbance and urbanization, were the 
primary factors contributing to its endangered status. On the one hand, it was recommended to 
strengthen the in-situ protection of the H. triflora, employing techniques such as artificial 
pollination to enhance gene flow among populations and thereby increase genetic diversity. On the 
other hand, we should focus on protecting wetlands from destruction. 

Key words: habitat fragmentation, RAD-seq, genetic diversity, genetic structure, species 
distribution models, demographic historical dynamics 
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保护 基因 组 学 (conservation genomics) 广义 上 是 指 将 基因 组 技术 和 方法 应 用 于 保护 生 
物 多 样 性 和 物种 保护 工作 中 的 一 门 学 科 CAllendorf et al., 2010) ， 主 要 聚焦 于 物种 的 濒危 机 
制 与 保护 策略 研究 .评估 物种 的 遗传 多 样 性 水 平 及 其 在 个 体 和 群体 中 的 分 布 格局 是 对 该 物种 
保护 遗传 研究 的 重要 内 容 (Hamrick & Godt, 1996) 。 分 子 遗传 标记 是 评价 种 质 资源 遗传 多 
样 性 的 重要 工具 ， 推 动 着 保护 遗传 学 的 发 展 。 虽 然 遗 传 信息 早已 运用 在 濒危 物种 的 研究 中 ， 但 
是 全 基因 组 范围 的 大 量 信息 标记 却 是 最 近 20 年 才 开 始 使 用 (Supple & Shapiro, 2018) ， 最 初 
也 只 是 应 用 在 少数 的 模式 植物 中 。 

高 通 量 测 序 技 术 Chigh throughput sequencing) 的 飞速 发 展 催生 了 简化 基因 组 测序 技术 。 
RAD-seq (restriction site Associated DNA) 是 一 种 常见 的 简化 基因 组 技术 ， 通 过 测序 限制 性 
beo 点 附近 的 短片 段 ， 具 有 和 覆盖 度 较 高 、 测 序 成 本 低 等 优点 。 它 不 受 参考 基因 组 的 限制 ， 

能 获得 大 量 的 高 质量 (single-nucleotide polymorphisms, SNP) 信息 位 点 (Davey et al., 2011; 
Narum et al., 2013) 。 目 前 RAD-seq 已 广泛 应 用 于 群体 研究 与 分 子 育 种 等 领域 ， 包 括 分子 标 
记 开 发 、 群 体 遗 传 分 析 、 遗 传 图 谱 构 建 和 全 基因 组 关联 分 析 等 。 孙 维 悦 等 (2022) 通过 简化 
基因 组 测序 对 濒危 植物 荷 叶 铁 线 蕨 (4diartura nelumboides) 的 保护 基因 组 学 研究 ， 揭 示 了 


其 濒危 机 制 。Cai 等 (20210 对 具有 重要 经 济 价值 的 热带 雨林 树种 演 南 风 吹 楠 CHorsfieldia 
tetratepala) 进行 RAD-seq 测序 ,制定 了 更 有 效 的 保护 策略 。 曹 航 鞭 等 (2022) 利用 RAD-seq 
对 极 和 危 物 种 波 叶 杜 鹏 (Rhododendron hemsleyanum) 的 遗传 特征 进行 无 参考 基因 组 分 析 ， 证 
明 对 于 无 参考 基因 组 的 极 小 种 群 野生 植物 (plant species with extremely small populations, 
PSESP) 也 可 用 简化 基因 组 测序 手段 获取 SNP 分 子 标记 。 因 此 ， 相 信和 针对 无 参考 基因 组 的 
濒危 物种 ， 简 化 基因 组 测序 可 以 获得 充足 的 遗传 信息 ， 为 保护 基因 组 学 研究 提供 有 效 工具 。 
此 外 ,有 效 种 群 大 小 的 估计 对 于 保护 和 管理 濒危 物种 、 了 解 种 群 动态 以 及 预测 遗传 演化 的 方 
向 和 速度 非常 重要 。 种 群 历史 动态 是 保护 基因 组 学 的 重要 研究 内 容 (Girod et al., 2011) 。 在 
种 群 的 演化 过 程 中 , 气候 变化 和 地 质 事件 作为 重要 的 外 部 作用 力 对 种 群 的 动态 变化 起 着 关键 
作用 


除了 物种 的 自身 特征 、 气 候 历 史 等 因素 影响 着 遗传 多 样 性 以 外 ,生境 片段 化 也 是 造成 生 
物 多 样 性 丧失 的 重要 因素 之 一 (Wilcove et al.,1998) ， 会 导致 遗传 多 样 性 降低 、 近 交 自 交 衰 
退 、 基 因 流 减少 和 种 群 规模 减少 等 (Young et ab, 1996; 陈 小 勇 ，2000; Honnay & Jacquemyn, 
2007) 。 还 会 造成 物种 栖息 地 减少 ， 影 响 传粉 者 数量 ， 降 低 传 粉 率 和 繁殖 成 功率 ， 从 而 影响 
g 物种 的 多 样 性 (Aguilar et al., 2006; González-Varo et al., 2009) 。 生 境 片 段 化 是 指 原先 连续 
> 的 生境 被 分 割 成 不 同 规模 和 形状 的 片段 而 形成 的 新 的 生境 格局 (Wilcove & Mclellan, 1986) 。 
目前 , 大 部 分 的 研究 表明 生境 片段 化 对 遗传 多 样 性 产生 的 负 效 应 多 于 正 效应 ( 武 唱 和 刘 志 民 ， 
2014) 。 在 以 专 性 传粉 的 石竹 (Dianthus deltoides) 研究 中 ， 发 现 片段 化 生境 的 石竹 其 种 子 
结实 率 低 并 可 能 产生 自 交 种 子 而 进一步 降低 其 种 群 的 遗传 多 样 性 〈Jennersten, 1988) 。 濒 危 
植物 Bankia goodii 受 生境 片段 化 影响 , 其 路 边 的 物种 个 体 繁殖 力 降 低 (Lamont et al., 1993 ) 。 
因此 ， 当 濒危 物种 遭受 生境 片段 化 影响 时 ， 应 进一步 加 强 研究 ， 以 更 好 地 理解 生境 片段 化 如 
可 影响 植物 遗传 多 样 性 ， 揭 示 其 濒危 机 制 并 提出 有 效 的 保护 策略 。 
水 角 CHydrocera triflora) 系 凤 仙 花 科 (Balsaminaceae) 水 角 属 多 年 生 水 生 草 本 植物 ， 
产 中 国 海南 、 印 度 、 斯 里 兰 卡 、 泰 国 、 越 南 、 老 挝 、 柬 埔 寨 、 马 来 西亚 及 印度 尼 西 亚 。 喜 生 
于 湖 边 、 沼 泽 湿 地 或 水 稻田 中 。 对 海南 省 水 角 资 源 现状 的 调查 发 现 ， 随 着 农村 旅游 的 开发 建 
设 ， 水 体 改 造 ， 湿 地 被 大 量 的 开发 利用 ， 水 角 生 境 受 到 不 同 程度 的 破坏 ， 该 种 分 布 范围 越 来 
越 狭 窄 , 十 分 罕见 ( 王 景 飞 等 , 2017) 。 在 IUCN 评估 报告 中 , 水 角 处 于 地 区 性 绝 灭 状态 CRE) 
级 别 。 有 关 水 角 的 相关 研究 较 少 ，2014 年 中 国 科学 院 昆明 植物 研究 所 在 海口 地 区 重新 发 现 
了 该 种 植物 ， 这 是 30 多 年 来 该 物种 在 国内 再 次 被 发 现 。 水 角 和 凤 仙 花 的 花 发 育 揭示 了 凤 仙 
花 科 最 早期 分 化 谱系 的 进化 趋势 , 花序 从 简单 的 轴 生 花序 逐渐 转变 成 更 复杂 的 腋生 花序 以 及 
花 被 片 的 分 化 是 一 个 逐渐 的 过 程 (Janssens et al., 2012) 。 水 角 作 为 凤 仙 花 科 的 单 种 属 ， 研 
究 其 遗传 多 样 性 不 仅 可 以 评估 生境 片段 化 对 遗传 多 样 性 的 影响 , 同时 也 为 凤 仙 花 科 的 进化 和 
系统 分 类 提供 基础 数据 。 
本 研究 以 来 自 海口 地 区 7 个 种 群 的 34 个 水 角 个 体 为 试验 材料 ,通过 RAD-seq 简化 基因 
组 测序 技术 ， 进 行 群体 遗传 分 析 ， 以 探 明 以 下 内 容 : CDD 基因 组 水 平 下 水 角 的 遗传 多 样 性 
和 遗传 结构 : OO 水 角 有 效 种 群 大 小 的 历史 变化 (3) 不 同 气候 情景 下 水 角 的 潜在 分 布 区 
变化 。 从 而 为 水 角 种 质 资 源 的 保护 利用 、 育 种 策略 的 实施 提供 理论 依据 。 


1 材料 和 方法 
1.1 实验 材料 

对 分 布 在 海口 市 的 水 角 进 行 野外 采集 (图 DD ， 共 收集 到 7 个 种 群 合 计 34 个 样本 。 
个 种 群 均 有 一 份 凭证 标本 ， 存 放 于 深圳 市 兰 科 植物 保护 研究 中 心 标本 馆 〈CNOCC) 。 每 个 种 
群 具体 信息 见 表 1。 
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表 1 水 角 样 品 采集 信息 
Table 1 Sample details of Hydrocera triflora 
ET Z5 BF Longitude 样本 量 
居 群 采样 地 点 纬度 M 
(E) : Individual 
Population Location Latitude (N) 
Number 
东兴 村 海南 省 澄 迈 东兴 村 Dongxing Village, 
110?10.7227 " 19?47.8291 " 5 
DXC Chengmai County, Hainan Province 
昌 旺 溪 海南 省 海口 遵 谭 镇 Zuntan Town, Haikou 
110°17.0243 " 19°47.2118" 5 
CWC City, Hainan Province 
卜 成 村 海南 省 海口 市 新 坡 镇 卜 茂 村 Bumao Village, 
110?18.7138 19?46.4791 " 4 
BMC Xinpo Town, Haikou City, Hainan Province 
EH 海南 省 海口 市 新 坡 镇 卜 壁 村 Bubi 
so Village, Xinpo Town, Haikou City, Hainan 110°20.7326” 19°46.8353 " 5 
Province 
bet 海南 省 海口 市 美 仁 坡 镇 卜 史 村 Bushi 
Age Village, Meirenpo Town, Haikou City, 110°23.4264" 19°49.6518" 5 
Hainan Province 
— 海南 省 海口 市 新 坡 镇 文选 村 Wenxuan 
Village, Xinpo Town, Haikou City, Hainan 110°20.0377” 19°44.7409” 5 
WXCI i 
Province 
Jobs 海南 省 海口 市 新 坡 镇 文选 村 Wenxuan 
M Village, Xinpo Town, Haikou City, Hainan 110?19.7120" 19?45.4986 " 5 


Province 


A. db: B. RSE; C. 生境 。 
A. Flowers; B. Fruit; C. Habitat. 


图 1 水 角 的 形态 特征 
Fig.1 Morphological characteristics of Hydrocera triflora species. 


1.2 DNA 提取 和 RAD 文库 构建 


本 研究 利用 天 根植 物 基因 组 DNA 提取 试剂 盒 CDP3050 进行 适当 改良 后 对 所 有 样品 的 
基因 组 DNA 进行 提取 。 通 过 琼脂 糖 凝 胶 电泳 和 Qubit 浓度 检测 后 ， 送 至 上 海 元 苹 生 物 医 药 
科技 有 限 公 司 。 采 用 EcoRI 限 制 性 内 切 酶 进行 酶 切 ， 然 后 构建 长 度 范围 在 300—500 bp 的 双 
端 测 序 文 库 。 纯 化 后 的 文库 在 Ilumina Novoseq 6000 测序 平台 完成 该 物种 简化 基因 组 测序 。 
1.3 获取 单 核 苷 酸 多 态 性 (SNP) 

在 测序 数据 下 机 后 ， 利 用 FastQC 软件 进行 质 控 ， 包 括 去 除 接头 、poly_A 和 低 质量 的 数 
据 ， 最 后 得 到 高 质量 的 数据 (Clean data) ， 用 于 后 续 的 分 析 。 本 研究 利用 Stacks v2.54 软件 

CCatchen et al., 2013) 对 水 角 的 RAD-seq 数据 进行 无 参 De Novo 组 装 和 变异 检测 。 分 别 运 

ÍT Stacks 软件 内 的 process. radtags 对 样本 进行 过 滤 , 除去 不 含 RAD 标签 的 序列 ,使 用 ustacks 
对 每 个 样本 建立 位 点 , 设置 控制 错 配 数 的 参数 M 为 3, 每 个 等 位 基因 所 需 的 最 小 片段 数量 m 
为 3; 运行 cstacks 程序 构建 包含 所 有 位 点 信息 的 catalog 目录 文件 ， 设 置 位 点 间 人 允许 的 错 配 
Zi n 为 3; 运行 sstacks 程序 把 每 个 样本 比 对 到 创建 的 目录 上 ; 再 运行 tsv2bam、gstacks 以 及 
population 程序 ， 得 到 包含 75 209 个 SNPs 位 点 的 vef 文件 。 为 保证 下 游 分 析 的 可 靠 性 ， 再 
使 用 VCFtools v0.1.13 软件 CDanecek et al., 2011) 对 Stacks 软件 中 检测 出 的 SNPs 进行 进 
步 的 过 滤 。 设 置 --max-missing 0.9， 即 最 大 缺失 率 为 90%; --maf 0.05， 将 最 小 等 位 基因 频率 
WEA 0.05; --minDPS， 确 保 基因 型 最 小 深度 为 5。 过 滤 后 剩 下 高 质量 的 38 667 个 SNPs 
用 于 后 续 群 体 遗传 结构 的 分 析 。 
1.4 群体 结构 分 析 
基于 34 个 水 角 样品 的 SNPs 位 点 , 运行 vef2phylip 脚本 构建 个 体 间 的 遗传 距离 矩阵 , 再 
利用 IQ TREE v1.6.9 软件 (Nguyen et al., 2015) 使 用 最 大 似 然 法 (maximum likelihood, ML) 
进行 系统 树 的 构建 。 根 据 1 000 次 重复 的 快速 自 展 (Ultrafast bootstrapping, UFBoot) (Minh 
et al., 2013) 和 SH-aLRT (Guindon et al., 2010) 测试 进化 校 ， 得 到 最 终 的 系统 发 育 树 。 使 用 
CN FigTree v1.4.3 软件 Chttp://tree.bio.ed.ac.uk/software/figtree/) 可视化 和 编辑 生成 的 树 文件 。 
c2 先 将 之 前 获得 的 高 质量 的 SNP 数据 集 使 用 VCFtools 软件 转换 成 ped 和 map 文件 , 再 将 
= 步 所 得 文件 进行 过 滤 生 成 bed 文件 。 最 后 用 Admixture (Alexander et al., 2009) 基于 每 个 


YT 


— 个 体 的 基因 型 推断 个 体 祖 先 成 分 , BEI A K 为 2 到 7 情况 下 的 分 群情 况 , 根据 不 同 K 
a 值 模拟 下 的 最 小 交叉 验证 〈cross validation, CV) 获得 最 佳 分 群 数 。 
本 研究 使 用 GCTA v1.26.0 软件 (https:/yanglab.westlake.edu.cn/software/gcta/) 生成 用 于 
PCA 分 析 的 和 矩阵， 再 用 了 及 中 的 ggplot2 软件 包 绘制 主 成 分 分 析 的 结果 图 。 
1.5 群体 遗传 多 样 性 分 析 
通过 Stacks 软件 中 的 Populations 子 程序 计算 物种 水 平 的 近 交 系数 Cinbreeding coefficient, 
Ps)， 观 测 杂 合 度 Cobserved heterozygosity, Ho) ， 期 望 杂 合 度 (expected heterozygosity, He). 
同时 基于 所 有 【变异 和 非 变 异 ) 位 点 计算 核 音 酸 多 态 性 (nucleotide diversity, Pi) 等 遗传 参 
数 。 再 用 R 语言 的 GroupGenome v2.6.1 程序 包 (Pfeifer et al, 2014) 计算 种 间 遗 传 分 化 系数 
(genetic differentiation, Fst) ， 衡 量 群 体 间 遗 传 分 化 水 平 。 
1.6 Mantel 相关 性 检验 和 分 子 方差 分 析 
为 了 评估 7 个 水 角 种 群 的 遗传 距离 、 地 理 距离 以 及 环境 距离 之 间 的 相关 性 ， 利 用 R 软 
件 中 的 “raster” 和 “rgdal” 提 取 每 个 种 群 点 的 气候 信息 ， 构 建 19 个 环境 因子 的 环境 距离 矩 
阵 ， 利 用 “geosphere” 构 建 地 理 距离 的 矩阵 ， 利 用 “ape” 读 取 遗 传 距离 矩阵 ， 再 利用 “ade4” 
对 两 两 矩阵 进行 9999 次 的 Mantel 重复 检验 , 最 后 用 “ggplot” 进 行 Mantel 检验 结果 的 绘图 。 
此 外 ， 利 用 Arlequin 3.5 软件 CExcoffier & Lischer, 2010) 进行 分 子 方差 分 析 (analysis of 
molecular variance, AMOVA) 。AMOVA 基于 进化 距离 〈evolutionary distance) 来 度量 并 计 
算 单 倍 型 〈 含 等 位 基因 ) 或 基因 型 间 的 方差 (820 ， 从 而 推断 群体 不 同 遗 传 组 分 之 间 的 变异 
水 平 。 


1.7 有 


效 群体 大 小 分 析 


ZE J 
于 无 参考 基因 组 的 物种 (Li & Durbin, 2011) 。 而 基于 SES 的 方法 不 依赖 于 位 点 之 间 的 连锁 
或 连续 性 CGutenkunst et al., 2009) ， 
Fu, 2015; Lapierre et al., 2017) 。 因 此 ，; 
CLiu & Fu, 2020) 1 
化 。 生 成 无 缺失 且 无 连锁 不 
用 python fi 
ID) SFS. 


SMC 模拟 种 条 


历史 动态 的 方法 是 高 度 依赖 于 


日 对 近期 的 


EW) 


IK ff 


将 SFS 状态 指定 为 folded， 以 统计 次 等 位 基因 


' 群 历史 大 小 变化 
各 所 获得 的 SNP 数据 使 


完整 基因 组 上 的 重组 图 谱 


的 SNP 数据 集 ， 将 VCF 文件 


， 并 不 适用 


TESTI HESS (Liu & 
基于 SFS 的 Stairway Plot v2 
fr 100 个 千年 前 (thousand years ago, kya) 到 1 kya 间 的 Ne 变 
F 衡 Clinkage disequilibrium, LD) 
JAX easySFS (https://github.com/isaacovercast/easySFS ) 建立 一 维 (one-dimensional, 


的 SFS。 选 择 使 位 点 数 尽量 多 的 


projection 值 输出 SFS 信息 ， 将 SFS 信息 输入 Stairway Plot v2 运行 所 需 的 blueprint 文件 中 。 


突变 率 设 置 为 1.0X103/ 位 点 / 代 ，2 年 为 一 代 ，200 次 


信 区 间 。 结 果 在 R 中 可 视 化 。 


1.8 18 


在 分 布 区 预测 


展 检验 用 于 计算 Ne 中 位 数 及 95% 置 


根据 材料 采集 ， 分 布点 经 纬度 信息 由 GPS 实地 定位 获得 和 通过 中 国 数字 植物 标本 馆 


Chttp://www.cvh.org.cn/) 和 全 球 生 物 多 样 性 
获取 标本 采集 地 经 纬度 信息 ， 共 收集 到 东南 亚 水 角 


zi ji Chttps://www.gbif.org/) 等 相关 数据 库 ， 


坐标 信息 200 余 个 。 对 所 有 数据 按照 


MaxEnt 模型 (http:/biodiversityinformatics.amnh.org/open_source/maxent) 要求 进行 标准 化 处 


理 ， 


并 且 为 了 防止 在 进行 最 大 烂 模 型 分 析 时 出 现 过 拟 合 现 象 ， 将 水 角 


的 分 布 数据 进行 得 选 ， 


在 50 km? 内 的 分 布点 仅 保留 一 个 , 共计 得 到 29 个 具有 经 纬度 信息 的 水 角 分 布点 并 整理 生成 


CSV 文件 。 


从 全 球 气 候 数据 库 WorldClim2.1 Chttp://www.worldclim.org) 分 别 下 载 全 新 世 中 期 MID, 
mid-Holocene; ~6000years BP) ~ “44% (1970—2000 F) 和 未 来 (2081—2100 Æ) 的 19 个 


生物 气 


(1.0 ©) 


适 生 


IREE, 


EEr 


J EE 


。 加 载 到 


2 结果 与 分 析 


2.1 测 
本 


FR fq) 格 式 储 存 , 包含 了 测序 的 序列 信息 和 其 对 应 的 测序 质量 信 
共 获 得 37.91 Gb 的 clean data 24 


919 到 


序数 据 概况 


4 


区 。 采 用 ROC I 


研究 利 


] RAD 简化 基因 组 技术 对 34 个 水 角 


55 093 740, 


分 布 的 主导 环境 因子 。 


数据 空间 分 辨 率 为 2.5 min (5k 分 辨 率 ) 。 未 来 气候 情景 选择 2 种 经 典 的 二 
氧化 碳 浓度 路 径 ， 分 别 为 RCP2.6 和 RCP8.5， 代 表 未 来 全 球 平 均 温 度 较 现在 增长 的 最 低 值 
直 (2.0 C). 
运行 MaxEnt 软件 模拟 水 角 在 不 
试 集 为 25% 进 行 模拟 分 析 。 将 物 
法 CGackknife) 模型 迭代 500 1X, itf 
并 制作 预测 图 
suitability index, HSI) ， 划 分 为 四 个 等 级 。 其 
区 ; 0.3 < HSI € 0.5 PEA 
结果 的 精准 度 ， 环 境 变量 的 贡献 率 判断 影响 物 


同 气 候 情 景 下 的 潜在 分 布 格局 ， 设 置 训练 集 为 75%， 测 
! 分 布 数据 与 气候 环境 数据 一 起 导入 MaxEnt 中 ， 选 择 刀 切 
各 环境 变量 对 物种 分 布 格局 的 贡献 率 
ArcGIS 10.2 中 ， 提 取 东 南亚 区 域 ， 根 据 各 
中 HSI < 0.1 为 非 适 生 区 ; 
EX; HSI > 0.5 为 高 适 生 


， 绘 制 响应 曲线 
E 境 适宜 度 指数 Chabitat 
0.1 < HSI x 0.3 为 低 
1 线 来 评价 模型 预测 


样本 进行 测序 ， 测 序 结果 以 FASTQ〔 简 
息 。 经 过 质量 评估 和 过 滤 后 ， 
E. 平均 每 个 样本 为 1.12 Gb. 每 个 样本 的 reads 读 长 在 2 512 
其 中 GC 含量 均值 为 36.2%， 所 有 样本 的 Q30 均值 为 92.9%， 表 明 测 序 


质量 较 高 ， 达 到 后 续 分 析 的 要 求 。 经 过 使 用 Stacks 进行 无 参考 基因 组 分 析 后 ， 得 到 水 角 34 
个 样本 的 平均 位 点 数 为 68 226， 平 均 片 段 覆 盖 深 度 为 8.61X。 经 过 VCFtools 再 次 过 滤 后 ， 

最 终 获 得 38 667 个 高 密度 的 SNP 位 点 。 
2.2 群体 遗传 结果 分 析 


进行 Structure 分 析 发 现 ， 


M2, 


K=3 It, CV 值 最 小 (图 2: A) ， 说 明 34 个 个 体 可 分 为 3 
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个 遗传 聚 类 (图 2: BO 。 利 用 最 大 似 然 法 (ML) 基于 38 667 个 SNPs 位 点 对 34 个 样本 构 
建 的 系统 进化 树 表 明 (图 2: CO , 每 个 种 群 的 个 体 都 能 聚 在 同一 个 遗传 分 支 上 , 并 具有 100% 
的 支持 率 。 昌 旺 溪 (CWC) . BEM (BMC) 、 东 兴 村 (DXC) 的 遗传 结构 更 相似 为 一 组 
(Clusterl) ; ÞE} (BBC) . bE (BSC)、 文 选 村 2 CWXC2) 的 遗传 结构 显示 出 更 
近 的 亲缘 关系 ， 倾 向 于 男 一 组 (Cluster2) ; 文选 村 1 OWXCI 表现 出 与 其 他 种 群 不 一 样 的 
遗传 结构 ， 为 Cluster3 。 主 成 分 分 析 (PCA) 是 根据 群体 的 原始 数据 矩阵 计算 出 单个 个 体 的 
特征 向 量 ， 再 绘制 特征 向 量 的 散 点 图 来 获得 群体 分 类 情况 。34 个 样品 基于 主 成 分 1 MER 
分 2 绘制 PCA 图 (图 3) 。 可 见 ， 每 个 种 群 都 分 布 得 很 集中 ， 这 与 系统 树 结果 一 致 。 其 中 
DXC, CWC, BMC, BBC 和 WXC2 分 布 相对 集中 。 在 PC2 维度 上 ，WXC1 与 其 他 种 群 相 
隔 得 较 远 ， 这 与 structure 分 析 的 结果 一 致 。 
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A. 不同 K 值 所 对 应 的 交叉 验证 错误 值 ; B. 最 佳 聚 类 结果 的 地 理 分 布 , 其 中 饼 状 图 代表 每 个 种 
AEA Tei; C. 基于 38 667 个 SNP 的 最 大 似 然 法 建树 结果 和 当 K-3 时 的 遗传 结构 分 析 结 果 。 
CH K-3 时 ， 每 一 个 柱 形 代 表 1 个 样本 ， 每 种 颜色 代表 1 个 遗传 秘 。) 


A. The admixture validation error rate corresponding to the different K values; B.The distribution 


of best results in clusters by Structure. The pie chart represents the genetic cluster of each 
population; C. Maximum likelihood tree and structural analysis based on 38 667 SNPs and Results 
of the genetic structure analysis! with K =3. (When K=3, each color represents a genetic cluster.) 


图 2 水 角 的 群体 遗传 结构 分 析 结果 


Fig. 2 Analysis of population structure 


颜色 分 别 代表 7 个 不 同 种 群 。 


Colors represent 7 different populations. 


population 
*- BBC 

n = BMC 

-= BSC 


PC2 


= CWC 

+ DXC 
Li » WXC1 
*- WXC2 


图 3 基于 SNPs 的 水 角 PCA 分 析 结 果 。 


Fig. 3 Results of the principal component analysis of Hydrocera triflora based on SNPs 


2.3 群体 遗传 多 样 性 分 析 


对 7 个 种 群 的 遗传 多 样 性 信息 进行 的 统计 分 析 表 明 〈 表 2) ，34 份 材料 的 3 


F 均 观察 杂 合 


E (Ho) 、 期 望 杂 合 度 (He) 和 核 昔 酸 多 样 性 〈 疡 ) 分 别 为 0.015 69. 0.165 4 和 0.186 5, 
因 最 高 (1445) ， 


同时 其 观察 杂 合 度 最 低 (0.079 1) 。 核 蔡 酸 多 样 性 (Pi) 分 析 表 明 ， 
Pi=0.01502~0.2399， 在 采样 地 东部 的 Cluster] 种 群 (BMC. CWC 和 DXC) 具 


各 种 群 中 均 具 有 特有 等 位 基因 , 平均 值 为 828.29， 其 中 BSC 的 特有 等 位 基 


核酸 多 样 性 


有 更 高 的 遗传 


多 样 性 。 近 交 系 数 Fis =-0.006 1~0.168 3， 种 群 内 近 交 系数 (Fis) 平均 值 为 0.061 67， 其 中 


CWC 和 WXC2 略为 负 值 ， 说 明 该 两 个 种 群 存在 杂 合 子 过 剩 现象 。 


基于 38 667 个 SNPs 位 点 的 34 个 水 角 样 本 种 群 间 成 对 的 遗传 分 化 系数 Fst 值 计算 结果 
( 表 3) 表明 : 遗传 分 化 系数 Fst 为 0.203 5-0.584 6， 其 中 CWC 与 BMC 之 间 的 Fsr 值 最 小 
(0.203 5) ， 表 明 两 个 种 群 间 的 遗传 交流 在 所 有 种 群 间 最 为 频繁 。 其 他 成 对 Fsr 均 大 于 0.25 


表明 水 角 7 个 种 群 间 的 遗传 分 化 程度 较 高 。 其 中 ， 文 选 村 1 (WXC1) 种 群 分 化 最 大 ， 与 其 


他 6 个 种 群 间 的 遗传 分 化 系数 均值 为 0.528 1。 
X2 水 角 7 个 种 群 的 遗传 多 样 性 


Table 2 Summary statistics of genetic diversity for seven populations 


种 群 名 称 ”特有 等 位 基因 ”观察 杂 合 度 。 EERE ” 近 交 系数 
Population Ap Ho He Pi Fis 
BMC 162 0.188 8 0.207 4 0.239 9 0.1 
CWC 389 0.210 5 0.184 1 0.206 5 -0.004 8 
DXC 1 385 0.179 8 0.177 0.199 3 0.041 3 
BBC 408 0.126 8 0.160 9 0.180 5 0.115 5 
BSC 1 445 0.079 1 0.143 8 0.161 0.168 3 
WXC2 756 0.172 3 0.150 1 0.167 8 -0.006 1 
WXCI 1 253 0.140 8 0.134 3 0.150 2 0.017 5 


平均 值 


828.29 0.156 9 0.165 4 0.186 5 0.061 67 
Average 
表 3 种 群 间 的 遗传 分 化 系数 
Table 3 Pairwise Fsr values between each population 
CWC BMC BBC BSC WXCI WXC2 
DXC 0.445 4 0.353 8 0.445 0 0.498 4 0.550 9 0.549 6 
CWC 0.203 5 0.404 9 0.497 4 0.539 6 0.520 6 
BMC 0.288 7 0.394 4 0.450 0 0.414 2 
BBC 0.392 1 0.498 7 0.5012 
BSC 0.544 9 0.548 7 
WXCI 0.584 6 


2.4 AMOVA 分 子 方 差分 析 和 Mantel 检测 


AMOVA 分 子 方差 分 析 的 结果 《〈 表 4) 显示 ， 当 把 7 个 水 角 
9.35% 的 遗传 变异 发 生 在 组 间 ， 群 组 间 的 遗传 分 化 为 Fer-0.093 53, THAT 


= 


不 显著 。 种 群 内 的 遗传 变异 占 总 变异 的 52.65%, 


离 和 环境 距离 的 Mantel 检测 ， 结 果 表 明 两 者 之 间 存 在 
(图 4: AO ; 基于 遗传 距离 和 地 理 距离 的 Mantel 检测 


J 


J 


群 分 成 上 述 3 个 群 组 时 ， 

f 组 间 的 遗传 分 化 
' 群 间 的 遗传 变异 率 为 38%。 基 于 遗传 距 
显著 正 相 关 (R=0.434 57，P=0.044 7) 


A (P=0.3393) (图 4: BO ; 基于 环境 距离 与 地 理 距离 的 Ma 


离 与 地 理 距离 存在 显著 的 正 相 关 关系 , 意味 着 两 者 之 


22 (图 4: ©). 


表 4 水 角 各 组 间 分 子 方差 分 析 


, d 


E 明 两 者 之 间 不 存 寿 
间 存 在 很 强 的 联系 (R=0.643 2, P=0.008 


Table 4 AMOVA analysis among the groups of Hydrocera trifloral. 


EXP 3 A ER 


自由 度 变异 组 分 变异 百分比 
变异 来 源 平方 和 | 
Degree of Variance Percentage of 
Source of variation Sum of squares 
freedom components variation 
群 组 间 
2 37 703.709 283.67204 Va 9.35 
Among populations 
组 内 种 群 间 
4 50 840.824 1152.49490 Vb 38 
Among populations within groups 
种 群 内 
61 97 408.6 1596.86230 Vc 52.65 
Within populations 
总 计 
67 185 953.132 3033.02923 
Total 
0.09353 


Fer 
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A， 成 对 遗传 距离 与 环境 距离 ，B， 成 对 遗传 距离 与 地 理 距离 ; C， 环 境 距离 与 地 理 距离 。 
A. Pairwise genetic distance (Fsr/1-Fsr) is associated with environmental distance; B. Pairwise 


genetic distance is associated with geographical distance; C. Environmental distance is associated 


with geographical distance. 
4 遗传 距离 、 地 理 距 离 和 环境 距离 之 间 的 Mantel 检测 结果 
Fig. 4 Mantel test results among genetic distance, geographi distance and environmental distance 
2.5 群体 大 小 历史 分 析 
为 在 Stairway Plot 中 对 水 角 的 群体 大 小 随时 间 变 化 进行 推测 , 首先 筛选 到 不 连锁 的 
SNP 数量 38 667 个 。Stairway Plot 推测 该 种 群 100 kya 至 1 kya 的 有 效 群体 大 小 变化 
5) 。 如 图 所 示 在 11.7 kya 的 末次 冰期 Cast glacial period, LGP) 以 后 , 种 群 有 收缩 趋势 。 
在 MID 时 期 ， 特 别 是 在 4 kya-2 kya， 水 角 有 效 种 群 大 小 出 现 明显 的 下 降 。 
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10+ 


Hydrocera triflora 


1 10 100 
Time (1k years ago) 


粗 线 为 中 位 数 ， 阴 影 部 分 为 95% 置 信 区 间 。 
Thick lines represent the median, and shaded areas represent the 95% confidence intervals. 
5 Stairway Plot 基于 SFS 的 群体 大 小 变化 分 析 。 
Fig.5 Changes in population size inferred by Stairway plot using folded site frequency spectrum 
(SFS). 


2.6 不 同情 景 下 的 潜在 分 布 格局 

根据 Maxent 模型 AUC 值 均 大 于 0.952， 表 明 该 模型 预测 物种 适宜 分 布 区 的 结果 比较 可 
靠 。 在 环境 变量 中 , 温度 季节 性 、 最 冷 月 份 的 最 低温 度 、 年 平均 温度 和 海拔 所 占 贡献 率 较 大 ， 
表明 它们 是 影响 水 角 分 布 的 主要 环境 因子 。 
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结果 表明 从 MID 时 期 到 当代 (图 6: A、6: B 及 表 5) ， 水 角 的 总 适 生 面积 大 幅度 的 
减少 , 主要 表现 在 高 适 生 面积 减少 29%, 中 适 生 面积 减少 34%。 未 来 在 二 氧化 碳 浓度 RCP2.6 
路 径 时 期 (图 6: CO ， 水 角 的 潜在 分 布 总 适 生 面积 较 当 代 增 加 了 4.74%， 主 要 表现 在 低 适 
生 面 积 较 当 代 增 加 了 33.40%。 在 二 氧化 碳 浓度 RCP5.8 路 径 时 期 (图 6: DO ， 水 角 的 总 适 
生 面 积 变动 不 大 ， 但 是 高 适 生 面积 大 幅 减 少 ， 高 适 生 面积 较 当 代 减 少 了 23.32%， 中 适 生 面 
积 较 现在 减少 了 4.91%， 而 低 适 生 面 积 增加 了 72.85%。 可 以 看 出 ， 在 未 来 ， 随 着 二 氧化 碳 
浓度 的 升 高 ， 全 球 平均 温度 的 升 高 ,海平 面 的 上 升 ， 生 境 破碎 化 严重 ， 对 水 角 的 高 适 生 区 有 


很 大 的 影响 ， 大 量 的 高 适 生 区 将 会 丧失 或 转化 为 低 适 生 区 。 


D. KÆ (2081—2100) RCPS. 5 气候 情景 。 


MID 


2081-2100 RCP2.6 
A. 全 新 世 中 期 情景 ; B， 当 代 (1970 一 2000) 气候 情景 ; C， 未 来 (2081 一 2100) RCP2. 6 气候 情景 ; 


1970-2000 


2081-2100 RCPS.8 


(图 中 虚线 圆圈 部 分 为 海南 岛 潜在 分 布 区 的 放大 图 。 深 


绿色 : 非 适 生 区 ; 浅 绿色 : 低 适 生 区 ; 黄色: PEEK; 红色: 高 适 生 区 ) 

A. mid-Holocene climate scenarios; B.Contemporary (1970 一 2000) climate scenarios; C.Future (2081 
—2100) RCP2.6 climate scenarios; D. Future (2081—2100) RCP8.5 climate scenarios. (The dashed 
circle in the diagram represents an enlarged view of the potential distribution area of Hainan Island. Dark 


green: Unsuitable; Light green: Low suitable; Yellow: Medium suitable; Red: Highly suitable) 
6 不 同情 景 下 水 角 的 潜在 分 布 格局 。 


Fig.6 Distribution of suitable habitats for Hydrocera triflora under different scenarios. 


X 5 水 角 在 不 同时 期 不 同 环境 变化 下 的 适 生 面积 


Table 5 The suitable area of Hydrocera triflora under different environmental changes in 


different periods 
总 适 生 面 高 适 生 中 适 生 低 适 生 不 适 生 
情景 积 Highly suitable Medium suitable Low suitable Unsuitable 
Climate Total 比例 
l , 面积 比例 面积 比例 面积 比例 面积 , 
scenarios suitable Proportio 
Area Proportion Area Proportion Area Proportion Area 
area n 
全 新 世 中 期 
590.61 116.84 3.07 205.76 5.4] 268.01 7.05 3 212.48 84.47 
MID 
现代 
416.75 82.95 2.18 135.65 3.57 218.91 5.76 3 365.80 88.50 


Contemporary 


未 来 RCP2.6 
Future RCP2.6 
未 来 RCP5.8 
Future RCP5.8 


436.51 86.55 2.28 154.05 4.05 292.02 7.68 3 270.69 86.00 


403.76 63.59 1.67 128.99 3.39 378.58 9:95 3 399.55 89.38 


3 讨论 与 结论 
3.1 水 角 的 遗传 多 样 性 : 生境 片段 化 降低 了 水 角 的 遗传 多 样 性 
自 20 世纪 90 年 代 以 来 ,生境 片段 化 问题 在 国内 逐渐 受到 关注 ， 一 直 是 保护 生物 学 和 生 
态 学 领域 的 研究 热点 之 一 。 起 初 集中 在 热带 雨林 中 开展 ( 朱 华 等 ，2001; ESÈ, 2016), 
亚热带 森林 次 之 〈 李 铭 红 等 ，2008) 。 在 不 同 生境 类 型 中 ， 岛 屿 的 研究 仅 占 2%《 吴 倩 倩 等 ， 
2017) 。 研 究 物 种 的 遗传 多 样 性 ， 不 但 可 以 反映 物种 适应 环境 变化 的 潜力 ， 而 且 还 可 以 为 物 
种 资源 特别 是 濒危 物种 现 有 的 生存 方式 及 状态 进行 有 效 的 评估 (Mable, 2019) 。 研 究 发 现 ， 
由 于 遗传 漂 变 和 近 交 的 原因 , 分布 广泛 的 物种 遗传 多 样 性 高 于 地 理 分 布 狭窄 的 稀有 和 濒危 物 
v Tit CCole, 2003) , 遗传 多 样 性 水 平 较 低 是 一 个 物种 濒危 或 分 布 狭窄 的 普遍 特征 (Gitzendanner 
= & Soltis, 2000; 景 昭 阳 等 , 2023) . AERE FE VE AL t Bt BE AE P L EE BE DISTR CCatchen 
et al., 2013) 。 本 研究 利用 SNP 分 子 标 记 对 水 角 7 个 种 群 进行 群体 遗传 学 分 析 ， 结 果 显 示 水 
具有 较 低 的 遗传 多 样 性 CPi-0.186 5; Ho-0.156 9; He-0.165 4) 。 在 濒危 植物 荷 叶 铁 线 其 
的 群体 遗传 分 析 中 也 有 类 似 的 发 现 ， 其 不 同居 和 群 的 平均 核 苷 酸 多 样 性 Pi、 观 察 杂 合 度 Ho 和 
期 望 杂 合 度 He 分 别 为 0.373、0.106 和 0.135《〈 孙 维 悦 等 ，2022) 。 与 木 本 植物 相 比 ， 草 本 
植物 具有 更 低 的 遗传 多 样 性 (Chung et al., 2020) 。 利 用 RAD-seq 技术 研究 的 极 危 树种 波 叶 
杜 鹏 的 遗传 多 样 性 为 (Pi=0.241 0; He=0.226 7) (HAS, 2022) 。 与 同 科 的 物种 相 比 ， 
水 角 的 平均 He 为 0.165 4 略 低 于 大 旗 瓣 凤 仙 花 Cmpatiens macrovexilla) 的 HeC0.018 6)(Zhou 
et al., 2020) 。 同 时 本 研究 结果 发 现 水 角 的 平均 观察 杂 合 度 小 于 平均 期 望 杂 合 度 ， 尽 管 结果 
在 统计 上 不 显著 , 但 这 也 暗示 着 水 角 种 群 中 存在 杂 合 子 缺 失 和 一 定 程 度 的 近 交 ,是 生境 片段 
化 的 遗传 效应 的 表现 。Dixo 等 (2009) 研究 发 现 ， 对 巴西 大 西洋 沿岸 森林 中 峰 肉 种 群 遗传 
s 多 样 性 与 破碎 化 面积 呈正 相关 ,生境 片段 化 对 剩余 栖息 地 之 间 的 基因 流 也 产生 负面 影响 。 生 
= 境 片 段 化 导致 种 群 之 间 的 遗传 流动 受到 限制 。 这 样 的 限制 可 能 导致 种 群 内 部 的 近亲 交配 增加 ， 
ET 杂 合 子 缺失 的 发 生 ， 进 而 降低 其 遗传 多 样 性 。 
对 不 同 种 群 间 的 遗传 距离 的 研究 有 助 于 对 不 同 种 群 的 遗传 多 样 性 进行 研究 , 也 可 用 来 表 
述 群 体 遗 传 结构 差异 。 对 于 不 同 种 群 之 间 的 群体 遗传 多 样 性 可 通过 遗传 分 化 系数 Cs. 来 
进行 研究 。 当 Fsr 为 0.25 以 上 ， 表 示 种 群 间 有 很 大 的 遗传 分 化 CWright, 1978) 。 在 水 角 7 
个 种 群 两 两 之 间 的 21 对 Fsr 中 有 20 个 Fst KF 0.25, 表明 7 个 种 群 间 的 遗传 分 化 水 平 较 高 。 
AMOVA 分 析 表 明 ， 水 角 38% 的 遗传 变异 发 生 在 组 内 群体 间 ， 而 群体 内 的 遗传 变异 更 高 ， 
为 S2.65%。 当 前 并 未 见 有 对 水 角 遗 传 多 样 性 的 研究 ， 凤 仙 花 科 仅 仅 包 括 水 角 属 和 凤 仙 花 属 
Qmpatiens) ， 这 些 植物 在 分 类 学 上 密切 相关 。 因 此 ， 对 凤 仙 花 科 的 研究 结果 可 以 为 水 角 的 
遗传 多 样 性 研究 提供 一 些 初步 的 参考 和 推测 。 钟 云 芳 等 (2014) 利用 ISSR 对 凤 仙 花 科 中 的 
研究 中 ， 发 现 海南 凤 仙 花 自 然 种 群 的 主要 变异 来 自 种 群 内 〈92%) ， 种 群 间 变 异 较 少 C890), 
水 角 也 有 与 之 一 致 的 结果 ， 即 种 群 内 的 遗传 变异 较 大 。 同 时 在 对 大 旗 办 凤 仙 花 的 研究 中 也 发 
现 ， 无 论 是 利用 ISSR 分 析 〈 李 祥 军 ，2013 ) 还 是 使 用 基于 “跳跃 基因 ”或 基于 逆转 录 转 座 
子 的 引物 间 结 合 位 点 (Inter-Primer Binding Site, iPBS) 标记 系统 来 评估 (Zhou et al., 2020) , 
其 种 群 内 的 遗传 多 样 性 大 于 种 群 间 的 , 即 主要 的 遗传 多 样 性 存在 于 种 群 内 ， 表 明 该 种 在 收集 
种 质 的 每 个 地 点 都 有 广泛 的 基因 型 。 在 聚 类 分 析 、 遗 传 结构 分 析 和 PCA 分 析 中 ， 我 们 发 现 


c 


本 研究 收集 到 的 海南 水 角 群 体 可 分 为 3 个 具有 不 同 遗 传 组 分 的 组 ,尽管 存在 这 些 不 同 的 遗传 
组 分 ， 但 其 组 间 的 遗传 分 化 较 低 (Fcr=0.09353) 。 考 虑 到 水 角 的 生境 特征 以 及 人 类 建筑 、 
农田 等 人 为 屏障 的 存在 ， 我 们 推测 亚 居 群 的 形成 可 能 是 由 于 生境 隔离 所 致 。 
对 水 角 的 环境 距离 、 遗 传 距离 和 地 理 距离 的 Mantel 相关 性 检验 结果 表明 : 水 角 环 境 虽 
离 与 地 理 距离 具有 显著 的 正 相 关 〈P=0.0082) ， 同 时 环境 距离 和 地 理 距离 也 存在 显著 的 正 相 
X (P=0.0447) ， 意 味 着 水 角 种 群 之 间 的 遗传 分 化 与 它们 之 间 的 环境 差异 有 关 ， 并 且 这 种 
环境 差异 与 地 理 距 离 的 变化 密切 相关 。 但 其 地 理 距 离 与 遗传 距离 没有 显著 相关 性 (P=0.339 
3) ， 即 水 角 种 群 的 遗传 距离 没有 随 着 隔离 距离 的 增加 而 增加 。 这 与 海南 凤 仙 花 自 然 种 群 的 
遗传 距离 和 地 理 距离 时 显著 的 高 度 相 关 性 (y=0.952，P<0.01) 有 差异 〈 钟 云 芳 等 ，2014) 。 
这 可 能 是 水 角 不 同 于 广泛 的 水 生 植 物 ， 其 具有 肉质 假 浆 果 , 裸露 于 水 面 或 沉 入 水 中 ， 可 供 鸟 
类 或 鱼 类 等 动物 取 食 , 因此 其 地 理 距 离 不 能 反应 其 种 群 的 遗传 距离 。 也 可 能 与 水 角 的 狭窄 分 
fü. 不 同 种群 间 其 地 理 距 离 过 近 有 关 。 总 的 来 说 ， 水 角 种 群 遗 传 分 化 水 平 高 ， 遗 传 距离 与 环 
境 距离 相关 但 与 地 理 距 离 不 相关 , 这 说 明 水 角 种 群 间 的 遗传 分 化 主要 受到 生境 片段 化 所 引起 
的 环境 距离 的 影响 ， 而 不 是 地 理 距离 的 影响 。 
3. 2 水 角 的 群体 历史 :万 年 前 琼 北 火 山 群 爆发 后 水 角 种 群 大 小 急剧 衰退 
~ 物种 的 种 群 历史 动态 是 其 所 经 历 的 历史 事件 和 遗传 特性 共同 作用 的 结果 (Hewitt, 2000), 
: 能 预测 物种 当前 和 未 来 的 生存 力 Hare et al., 2011) 。 本 研究 通过 使 用 大 量 的 SNP 标记 对 水 
角 进 行 群体 大 小 历史 分 析 , 发 现 水 角 在 100 kya 到 末次 冰期 之 前 , 有 效 种 群 大 小 都 相对 稳定 ， 
我 们 推测 这 与 全 新 世 早期 ， 全 球 性 气候 转 温 ， 海 南 岛 植 物 区 系 较为 兴盛 发 达 有 关 〈 颜 家 安 ， 
2006) 。 水 角 的 有 效 种 群 大 小 从 末次 冰期 (LGP: 11.2 kya) 后 开始 持续 下 降 ， 直 至 现在 〈 图 
5) 。 琼 北 火山 群 的 喷发 活动 可 能 发 生 在 距 今 约 1 AAA ERRE, 2004) ， 这 与 水 角 
有 效 种 群 大 小 开始 下 降 时 间 吻 合 , 推测 琼 北 火 山 群 的 喷发 活动 可 能 引起 土地 破坏 、 气 候 变化 
或 其 他 环境 因素 的 改变 ， 导 致 水 角 生 境 的 片段 化 ， 限 制 了 水 角 的 基因 交流 。 在 MID 时 期 ， 
水 角 有 效 种 群 大 小 在 明显 下 降 。 全 新 世 中 期 ， 是 一 个 全 球 范 围 内 人 类 活动 扩张 的 重要 时 期 。 
在 这 个 时 期 ， 人 类 社会 经 历 了 从 狩猎 采集 生 活 方式 向 农业 和 定居 生活 的 转变 , 这 对 环境 演化 
也 产生 了 重要 影响 。 在 新 石器 中 晚期 ， 有 0 万 余人 现 食 于 海南 岛 ， 与 植物 争 地 ， 植 物 生境 也 
逐渐 发 生 较 大 的 变化 《司徒 尚 纪 ，1987) . HAA (2002) 对 海南 岛 北 部 火山 口 沉积 孢 粉 
e 的 研究 ， 发 现 海南 岛 的 农业 活动 早 在 3000aBP 就 开始 对 地 低 植物 具有 较 明 显 的 影响 。 由 此 
一 推测 水 角 的 有 效 种 群 大 小 在 历史 时 期 首先 受到 琼 北 火山 群 爆发 的 影响 , 进而 又 受到 人 类 活动 
oO 的 影响 ,两 种 因素 又 加 导致 生境 的 持续 变化 。 这 种 有 效 种 群 大 小 的 减少 可 能 导致 了 其 遗传 多 
样 性 的 丧失 和 种 群 的 遗传 分 化 ， 逐 渐 形 成 如 今 的 分 布 格局 。 
3.3 水 角 的 潜在 分 布 格局 : 高 C0, 排 放 情 境 下 水 角 的 高 适 生 区 转变 为 低 适 生 区 
在 环境 因子 贡献 率 的 分 析 中 ， 本 研究 分 析 所 得 , 温度 季节 性 、 最 冷 月 份 的 最 低温 度 和 年 
平均 温度 的 总 贡献 率 占 比 达 67.6%， 这 与 水 角 喜 温 怕 寒 ， 温 度 低 于 15 CC 即 停止 生长 的 生活 
习性 相符 〈 王 景 飞 等 ，2017) 。 而 海拔 占 了 第 四 高 的 贡献 率 7.4%， 表 明 它 是 影响 水 角 分 布 
的 主要 环境 因子 之 一 , 也 佐证 了 在 海南 凤 仙 花 的 研究 中 低 海 拔 的 遗传 多 样 性 远 远 高 于 高 海拔 
梯度 的 种 群 这 一 结果 ， 说 明海 拔 是 影响 群体 基因 流 的 重要 因素 〈 钟 云 芳 等 ，2014) 。 
对 不 同 气候 情景 下 东南 亚 地 区 的 水 角 的 潜在 分 布 适 生 区 进行 模拟 分 析 , 结果 表明 从 MID 
时 期 到 当代 ， 水 角 的 总 适 生 面积 大 幅度 的 减少 C50 ， 其 中 位 于 海南 岛 的 适 生 区 由 高 适 生 
区 和 中 适 生 区 的 分 布 格局 转变 为 当代 的 以 低 适 生 区 为 主 的 分 布 格局 〈 图 6: A、6: B) X 
与 种 群 历史 大 小 的 分 析 结 果 一 致 ， 表 明 从 MID 时 期 开始 ， 水 角 的 生存 受到 威胁 直至 现在 。 
与 当代 气候 相 比 , 在 未 来 气候 变化 下 水 角 的 潜在 分 布 区 总 面积 变动 不 大 , 但 低 适 生 面积 较 当 
代 增 加 了 33.40%。 特别 是 在 二 氧化 碳 浓 度 RCP8.5 路 径 的 气候 环境 情景 下 ,高 适 生 面积 较 现 
在 减少 了 23.32%， 取 而 代 之 的 是 低 适 生 面积 增加 了 72.85%。 其 中 ， 位 于 马 来 群 岛 的 适 生 区 
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几乎 完全 消失 (图 6: DO . 根据 IPCC Cintergovernmental panel on climate change, IPCC) 发 

表 的 《第 六 次 评估 报告 》， 在 SSP5-8.5 (shared socioeconomic pathway, SSP) 的 气候 环境 下 ， 

二 氧化 碳 排放 量 到 2050 年 将 增加 一 倍 ， 到 2100 年 全 球 平均 气温 将 高 4.4 摄氏 度 。 我 们 推测 

这 可 能 与 气候 变 暧 条 件 下 亚热带 北 界 向 北 推移 有 很 大 的 关联 ， 未 来 CO; 浓度 升 高 将 引起 全 

球 变 暖 ， 动 植物 的 分 布 逐 渐 向 高 纬度 和 高 海拔 地 区 移动 〈Tingley, 2014) ， 从 而 导致 水 角 在 

马 来 群 岛 的 适 生 区 近乎 消失 。 这 些 分 析 结 果 表 明 , 未 来 的 气候 变化 将 进一步 影响 水 角 的 适 生 

区 , 特别 是 高 适 生 区 将 面临 丧失 和 转变 的 风险 。 这 说 明 水 角 种 群 面临 着 适 生 环境 变化 的 挑战 。 

人 类 开荒 、 伐 木 、 放 牧 等 活动 ， 以 及 交通 和 水 库 建 设 等 是 造成 全 球 的 自然 生境 和 生物 多 
样 性 都 面临 着 生境 丧失 和 片段 化 威胁 的 主要 原因 (Laurance et al., 2014; Haddad et al., 2015) 。 
城市 化 也 加 剧 着 栖息 地 的 破碎 化 ， 特 别 是 对 于 城市 中 残余 的 小 种 群 本 土 植物 更 是 如 此 。 

对 于 水 角 而 言 ， 较 低 的 遗传 多 样 性 水 平 ， 可 能 削弱 了 其 在 迁徙 过 程 中 的 适应 潜力 ， 导 至 
其 有 效 种 群 规模 降低 。 种 群 内 的 遗传 分 化 水 平 较 高 ， 说 明 其 自然 种 群 之 间 的 基因 流动 较 少 。 
环境 距离 与 遗传 距离 、 地 理 距离 之 间 都 存在 显著 的 正 相 关 关 系 , 表明 环境 因素 在 水 角 种 群 遗 
传 分 化 中 起 到 重要 的 作用 , 而 地 理 距离 则 作为 间接 因素 通过 影响 环境 差异 来 影响 遗传 分 化 的 
发 生 。 因 此 限制 水 角 种 群发 展 的 重要 原因 来 自 于 其 自身 较 低 的 遗传 多 样 性 和 生境 因素 , 具体 
表现 为 湿地 生境 的 片段 化 严重 , 再 加 上 人 为 活动 干扰 、 城 市 化 等 不 利 的 环境 条 件 使 得 其 生境 
进一步 被 破坏 ， 种群 规 模 进 一 步 减 小 ， 逐 步 形 成 如 今 地 区 性 绝 灭 状态 。 针 对 水 角 濒临 灭绝 这 
一 现象 ， 一 方面 要 着 重 保护 湿地 免 受 人 为 因素 的 破坏 ， 其 次 进行 就 地 保护 ， 更 有 效 的 提高 其 
群体 的 遗传 特性 。 另 一 方面 进行 人 工 授粉 等 方法 提高 其 基因 流 ， 以 促进 其 种 群 的 发 展 。 
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